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Abstract: The presence of a thick snowpack could interfere with forest 
stability, especially on steep slopes with potential damages for young and 
old stands. The study of snow gliding in forests is rather complex be¬ 
cause this phenomenon could be influenced not only by forest features, 
but also by snow/soil interface characteristics, site morphology, meteoro¬ 
logical conditions and snow physical properties. Our starting hypothesis 
is that different forest stands have an influence on the snowpack evolu¬ 
tion and on the temperature and moisture at the snow/soil interface, 
which subsequently could affect snow gliding processes and snow forces. 
The aim of this work is therefore to analyse the snowpack evolution and 
snow gliding movements under different forest covers, in order to deter¬ 
mine the snow forces acting on single trees. The study site is located in a 
sub alpine forest in Aosta Valley (NW-Italy) and includes two plots at the 
same altitude, inclination and aspect but with different tree composition: 
Larch (Larix decidua) and Spruce (Picea abies). The plots were equipped 
with moisture and temperature sensors placed at the snow/soil interface 
and glide shoes for continuous monitoring of snow gliding. The recorded 
data were related to periodically monitored snowpack and snow/soil in¬ 
terface properties. Data were collected during two winter seasons (2009- 
10 and 2010-11). The snow forces on trees were analytically calculated 
either from snowpack data and site morphology or also from measured 
snow gliding rates. Different snow accumulations were observed under 
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the two different forest stands, with a significant effect on temperature 
and moisture at the snow/soil interface. The highest snow gliding rates 
were observed under Larch and were related to rapid increases in mois¬ 
ture at the snow/soil interface. The calculated snow forces were generally 
lower than the threshold values reported for tree uprooting due to snow 
gliding, as confirmed by the absence of tree damages in the study areas. 

Keywords: snow/soil interface; temperature; water content; snow forces; 
trees 

Introduction 

Snow forces are the result of glide and creep processes of the 
snowpack: a) snow gliding is the slow downhill motion of the 
snow on the ground mainly influenced by the roughness of the 
ground surface and the wetness of the lowermost boundary layer 
of the snow cover; b) snow creep is the result of settlement and 
internal shear deformation parallel to the slope (In der Gand and 
Zupancic 1966). 

The snow gliding is strictly dependent on physical snow char¬ 
acteristics, overall depth, density and moisture, and on morpho¬ 
logical features, overall slope angle and surface roughness 
(McClung and Larsen 1989; Jones 2004; Margreth 2007; Holler 
et al. 2009). In the first period of extensive research of glide 
processes and glide avalanches, Bader at al. (1939) found that the 
most important factors that have an effect on glide are slope, ex¬ 
posure, degree of surface roughness, temperature at the ground 
surface and the thickness and properties of the snow cover. They 
proposed that the primary action of these factors is their effect on 
the frictional conditions of the slide surface. Basal friction is ob¬ 
viously reduced by the presence of liquid water at the snow/soil 
interface (McClung and Clarke 1987). Also Clarke and McClung 
(1999) stated that gliding can occur if 1) the interface is smooth 
(e.g. bare rock or grassy vegetation), 2) the temperature at the 
snow/ground interface is at 0° C, guaranteeing the presence of 
free water at the interface, and 3) the slope angle is at least 15° 
for roughness typical of alpine ground cover. Moreover, the in- 
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tensity of gliding rates is strictly related to increases in snow/soil 
interface temperature and moisture (Lackinger 1987; McClung 
and Larsen 1989; Clarke and McClung 1999; Newesely et al. 
2000 ). 

The wide-scale analysis of snow gliding intensity in alpine and 
subalpine areas has to take into account the land-use changing: 
the abandonment of pastures and the natural reforestation could 
cause a significant change on the ground roughness. In particular, 
Newesely et al. (2000) found that the abandonment of agriculture 
land may foster an increase of snow gliding rates due to the 
reduction of basal frictions. 

Although forests can reduce the snow gliding intensity 
(Leitinger et al. 2008) by preventing snow accumulation and 
homogeneous snow distribution due to interception and falling of 
snow from the crowns (Holler 1995, 2001), it is not possible to 
totally exclude snow gliding in forest. Snow gliding might have 
significant effects on forest plants; in particular, juvenescent 
trees might be uprooted from the ground or the growth of young 
plants might be strongly affected by the snow forces (In der 
Gand 1978). According to these evidences, Holler et al. (2009) 
evaluated, with different approaches, the snow gliding rates in an 
Austrian subalpine stand of small trees (Pinus cembra, Larix 
decidua and Picea abies). In the same work, snow forces 
measured in field, and back-calculated with different equations, 
were compared to the forces necessary to uproot the young trees: 
the authors reported that forces equal to 950 N and 3500 N can 
uproot juvenescent trees with a diameter of 0.02 m and 0.045 m, 
respectively. 

The study of the snow gliding phenomenon in forests is rather 
complex because it is influenced not only by forest features 
(species, crown, diameter and density), but also by topography 
(slope angle, morphology, aspect), meteorological characteristics 
(air temperature, rain on snow events) and snowpack structure. 
Moreover, the impact of snow forces on trees depends on the 
forest characteristics itself, and could lead to different 
consequences in term of tree stability. 

The aim of this work was to analyse the snowpack evolution 
and snow gliding movements under different forest covers, 
namely Larch (Larix decidua) and Spruce (Picea abies) stands, 
in order to determine the snow forces acting on trees of different 
species and diameters. In particular, we aim at evaluating the 
influence of physical properties as temperature and moisture at 
the snow/soil interface on snow gliding processes, and 
consequently on snow forces. 

Materials and method 

Study site 

As this research was developed within the project “Forets de pro¬ 
tection: techniques de gestion et innovation dans les Alpes occi¬ 
dentals” (Protection forests: management and innovation tech¬ 
niques in the western Alps) in the frame of the Operational Pro¬ 
gram ALCOTRA Italy-France (2009-2012), the study site was 
chosen within a snow avalanche protection forest. Located in the 
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Municipality of Brusson, in the Aosta Valley Region (North 
Western Italian Alps), at an elevation of 1950 m a.s.l., the site is 
characterized by a south aspect and a mean slope angle of about 
40°. The coniferous mixed forest of the study site is mainly 
composed by larch (Larix decidua ), spruce (Picea abies) with a 
lower presence of pine (Pinus uncinata). The lower part of the 
forest was reforested in 1940, in order to protect the village be¬ 
neath, and it was grazed. 

The climate of the area is continental humid subarctic, with a 
mean annual air temperature equal to ± 4 °C and a mean annual 
precipitation equal to 718 mm. The mean annual snow depth is 
175 cm (historical dataset from 1913 to 2001, Mercalli et al. 
2003). 

Two adjacent plots (40x40 m), located at the same elevation 
and aspect, were chosen (Fig. 1): a) one within a larch stand 
(Larch plot) and b) one within a spruce stand (Spruce plot). By 
counting all trees in each plot (with stem diameter > 5 cm) the 
average stem density was determined equal to 300 stems/ha for 
both plots. The crown canopy, measured as crown projection, 
amounted to 60% for the Larch plot and 65% for the Spruce plot. 
The average stem diameter was 40 cm and most of the trees were 
included in the range 30^10 cm (10 cm interval). The natural 
regeneration was generally scarce. 

Data collection 

At the snow/soil interface (2 cm deep in the soil organic horizon) 
each plot was equipped with thermistor probes, with accuracy of 
±0.1 °C, (MadgeTech®) and volumetric moisture probes, which 
measure the dielectric constant of the media through the 
utilization of capacitance/frequency domain technology, with 
accuracy of ±3% VWC, (EC-5MadgeTech). Temperature and 
moisture measured at snow/soil interface were named 
respectively as Ti and Mi, while the increase of Mi over time 
was named wetting rate. 

Moreover, two glide shoes connected to potentiometers 
(Sommer®) (Fig. 2) as developed by in In der Gand (1954) and 
successively applied by other scientists (e.g. Leitinger et al. 2008; 
Holler et al. 2009), were placed in each plot at a minimum 
distance of 5 m below the first upwards tree. The potentiometers 
were calibrated in order to convert the electric signal registered 
by the snow shoes into downward movement (mm/day), caused 
by snow gliding. All sensors were connected with a central data 
logger (Datataker® DT-80) supplied by a long duration battery 
and set to record the different parameters every 15 minutes. The 
electric wires were cabled underground in order to minimize 
their effect on the surface roughness. 

The meteorological data were acquired by an automatic snow 
and weather station (Gressoney-Weissmetten 2038 m a.s.l.) 
located 3 km apart from the study site (belonging to the Ufficio 
Centro Funzionale of the Aosta Valley Region). 

The snowpack physical characteristics (snow depth -H s ~, snow 
density -p-, grain types and size, snow temperatures, hand hard¬ 
ness index) were measured monthly, following standard interna¬ 
tional methods (Fierz et al. 2009): 8 snow profiles in winter 
2009-2010 and 3 snow profiles in winter 2010-2011 were de- 
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scribed in both plots (Table 1). The damages to the trees were and broken branches at the end of winters, 
evaluated in the field by observations of eventual uprooted trees 
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Fig. 1 Localization of the study site and of the two selected plots. 



Fig. 2 Moisture (a) and temperature (b) sensors, and a glide shoe (c), connected to the potentiometer, used in this study. 
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Table 1. Physical snow properties - according to Fierz et al. (2009) - measured in the field during winter 2009-2010 (8 snow profiles) and 
2010-2011 (3 snow profiles), under Larch and Spruce plots. 

Parameter 

2009- 

-2010 

2010- 

-2011 


Larch 

Spruce 

Larch 

Spruce 

Maximum value of snow depth (cm) 

60 

66 

28 

22 

Dominant grain type 

FCso, MFpc 

FCso, DHcp 

RGsr 

RGsr 

Maximum grain size (mm) 

2.5 

4 

1.5 

1.5 

Mean hand hardness Index 

3 

4 

3 

2 

Mean snow density (kg-m‘ ) 

277 

229 

220 

270 

Dominant grain type in the basal layer 

FCso, DHcp 

DHcp 

FCso 

FCso 

Mean snow/soil interface temperature (°C) 

-3.8 

-2.1 

-1.1 

-1.8 

Date of snowpack isothermal condition 

March 1 l rd 

March 19 rd 

March 14 rd 

March 20 rd 


Note: FCso means solid faceted particles; MFpc means rounded polycrystals; DHcp means hollow cups; RGsr means small rounded particles. 


Data collection was done from the middle of November 2009 
until the end of April 2010 (called winter 2009-2010) and from 
the middle of November 2010 until the end of April 2011 (called 
winter 2010-2011). The seasons were identified according to 
meteorological quarterly subdivision: fall included September, 
October, November; winter included December, January and 
February; spring included March, April and May. 

Calculation of snow forces 

The calculation of snow forces on stems was done following the 
approach used by Holler et al. (2009). In particular, two different 
methods were applied by using a) empirical equations (Margreth 
et al. 2007) and b) snow gliding data recorded with glide shoes 
(Me Clung and Larsen, 1989). We calculated snow forces for 
different stem diameters (0.02, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5 m) in 
order to consider the real diameter distribution observed in the 
field and to compare our results with the available literature. 

Snow pressure according to the Swiss guidelines 
Margreth et al. (2007) applied the Swiss Guidelines with some 
improvements in order to calculate snow pressure on cableway 
masts. The snow pressure was calculated with the following 
equation (Eq.l): 

S'n = pg(H s 2 / 2)KN (i) 

where: S N is the snow pressure component in the line of slope 
per meter run of the supporting surface along the contour line 
[kN-m 1 ], p is the density of the snow cover [t-m" 3 ], g is the 
gravitational acceleration [ms' 2 ], H s is the vertical snow depth 
[m]; K is the creep factor dependent on the slope inclination yj 
and the density p, and N is the glide factor. 

In Margreth et al. (2007) the K factor was associated to fixed 
snow density values. We interpolated those values in order to 
calculate the exact K factor from the snow density that we 
measured in the field (y = 0.73 x + 0.547, r 2 = 0.991, where y 
corresponds to K/srnly/ and x to p). After that, K/smlyj was 
multiplied by sin2 yj (where yj corresponds to the inclination of 
the study site). 
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In Margreth et al. (2007), the N factor depends on the ground 
roughness and aspect: it increases from WNW-N-ENE to ENE- 
S-WNW and it is inversely related to the ground roughness. 
Forest is not a determining parameters in the choice of N. 
Therefore, in our forested plots, we used the ground class 1, 
assuming that stems determined a high surface roughness (more 
than boulders). The study site was included in the ENE-S-WNW 
aspect class, because it is located on a South slope. Finally, we 
chose a glide factor N equal to 1.3. 

Equation 1, which is valid only for infinitely long planes, was 
integrated, as suggested by Margreth et al. (2007), with a 
coefficient ( rj F ), in order to consider the end-effect forces on 
narrow obstacles. q F was calculated with the following equation 
(Eq. 2): 

tjF = 1 + c(D/W) (2) 

where: D is the snow thickness (calculated as H s *cosy/, where H s 
is the vertical snow depth measured in the field and y/ is the slope 
angle); W is the width of the structure (that corresponded to the 
stem diameter); c is a non-dimensional coefficient and depends 
on the snow gliding intensity. It is usually deduced from the 
Swiss Guidelines. We used the smallest coefficient (c = 0.6) 
because the snow gliding is typically low in dense forest areas 
(Holler 2001). 

Subsequently, equations 1 and 2 were joined in Eq. 3: 

S'NM = pg(H s 2 / 2)KNti f (W / D) (3) 

> 

Finally, the snow pressure S N>M was multiplied by the different 
stem diameter classes in order to compare our results, in term of 
forces, with those presented by Holler et al. (2009). 

Snow pressure according to Me Clung and Larsen (1989) 

This method considered the snow gliding rates directly measured 
in the field to determine principal stress on structures in 
dependence of the stagnation depth: 
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<7 = {2 /(I- v)][(d 7 D) + (L/H S )\ }0.5/?gdisn f + 
0.5[ v /(I - v) /?gDcos y/\ 

where: o is the maximal principal stress of the snow on a 
structure [Nm 2 ], v is the Poisson coefficient deduced from Salm 
et al. (1971). We used v = 0.1 that corresponds to a mean snow 
density of 250 kg-m" 3 , equal to the density we measured in our 
plots; LIH S is a dimensionless creep parameter calculated by 
McClung and Larsen (1989) equal to 0.27 + (w/12), D is the 
snow thickness, p is the snow density measured in the field 

[kg-m' 3 ], g is the gravitational acceleration [m-s' 2 ], yj is the slope 

> 

angle, and d is the stagnation depth. 

The stagnation depth is a fundamental parameter to formulate 
the boundary condition at the snow/soil interface. We calculated 
the stagnation depth d’ from the measured snow gliding rates 
using Equation 5 (McClung 1974): 

d' = \jiv g ] 42(1 - v)pgD sin y/\ 

where: p is the viscosity, equal to 10 11 Pa for snow, V g is the 
snow gliding velocity [mm/day]. According to Hoeller (2009) we 
used the maximum daily snow gliding rate recorded during 
winters 2009-2010 and 2010-2011; v is the Poisson coefficient. 

Statistical analysis 

Descriptive statistics were used to analyzed meteorological 
pattern, snow physical characteristics, and snow gliding data. 
Correlation analysis between temperature and moisture at 
snow/soil interface and between air temperature and temperature 
and moisture at snow/soil interface were carried out using the 
SPSS 12.0 for Windows (SPSS, 2003). 

Results and discussion 

Winter 2009-2010 
Meteorological patterns 

The mean air temperatures was -5.2 °C with a minimum value of 
-14.3 °C registered on 19th of December 2009 (Fig. 3a). This 
season was characterized by one short mild period, from 
December 6th to 10th, when the air temperature reached positive 
values with a maximum of +3.9 °C. The mean daily air 
temperature became constantly positive after 20th of April 2010. 
The automatic weather station recorded a cumulative 
precipitation equal to 330 mm of snow water equivalent (SWE) 
for the whole winter. The maximum snow depth (135 cm) was 
reached on 2nd of April 2010. 

The cumulative snowfall was lower than the 75 years 
historical mean (488 mm of SWE) reported by Mercalli et al. 
(2003): in particular, the difference was remarkable especially 
during January, February and April. 


Snow physical characteristics 

In the Larch plot , the measured snow depth reached a maximum 
of 60 cm on January 13th.The snow density ranged between 184 
kg-m' (on January 13th) and 362 kg-m' (on March 3rd). The 
first snow profile (December 9th) showed an isothermal 
snowpack with a basal layer characterized by slush crystals 
(MFsl). From the middle of December to the middle of February 
the thermal gradient increased and determined the formation of 
faceted crystals (MCso) and depth hoar (DHcp). In March the 
snowpack was again mostly at isothermal conditions with a 
predominance of melt-freeze crusts (MFcr). 

The maximum snow depth measured in the Spruce plot (66 cm) 
was reached on February 2nd.The snow density ranged between 
163 kg-m' (on December 22nd) and 296 kg-m' (on March 3rd). 
From the middle of December till the middle of February the 
thermal gradient increased and determined the formation of 
faceted crystals (MCso) and depth hoar (DHcp). In March the 
snowpack reached isothermal conditions with a predominance of 
cluster rounded crystals (MFcl). 

Temperature and moisture at the snow/soil interface 
In the Larch plot , the temperature at the snow/soil interface 7) 
(Fig. 3b) was positive until December 10th, except on November 
29th, when a first freezing event occurred (daily average 7) = -1.3 
°C). Later 7} remained below 0 °C from the middle of December 
until the end of March. A sharp increase of 7) (which reached 0 
°C) was observed on December 26th, January 1st, February 19th, 
March 3rd and March 18th. The moisture at the snow/soil 
interface M, (Fig. 3c) was correlated with soil (r = +0.498; p 
<0.01) and air temperatures (r = +0.418;;? <0.01). The maximum 
values were recorded several times during this winter season and 
in particular on December 8th (21%), with a wetting rate of 
2%/day for 11 days, on December 25th (20%), with a wetting 
rate of 15% in one day and on February 26 th (22%) with a 
wetting rate of 5%/day for 2 days. The period between the 
beginning of January and the middle of February was 
characterized by a constant value of M, close to 10%. 

In the Spruce plot the temperature at the snow/soil interface 7) 
(Fig. 3b) was positive until November 29th (7)= -0.8 °C). Later 
f remained constantly below 0 °C from the middle of December 
until the end of March. The moisture at the snow/soil interface 
Mi (Fig. 3c) showed remarkable oscillations until the beginning 
of January, with 3 maximum values recorded on November 16 th 
(14%), with a wetting rate equal to 2.5%/day for 4 days, on 
December 11th (15%), with a wetting rate of 1%/day for 14 days 
and on December 25th (13%), with a wetting rate of 6.5%/day 
for 2 days. Later, M, remained constant (around 5%) until March 
19th. From March 19th to 24th a remarkable increase was 
observed (+3.2%/day), reaching a moisture value of 21.3% on 
March 24th. M, was correlated with 7) (r = +0.641; p <0.01) and 
air temperature (r = +0.631 ,p <0.01). 

Snow gliding 

While no snow gliding events were observed in the Spruce plot , 
under Larch , only a single event of snow gliding occurred (Fig. 
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3d) between December 1st and 14th after the first seasonal snow¬ 
fall (30 cm of fresh snow). In this period, Mi rose by 20% and 
the dominant grain type in the basal layer observed on December 
9 th was slush (MFsl), indicating a consistent presence of liquid 
water. The cumulative movement of the glide shoes was equal to 
19 mm, with an average daily snow gliding rate of 1.9 mm/day 
for 10 days and a maximum daily snow gliding rate of 3.9 


mm/day. The volumetric moisture at the snow/soil interface in¬ 
creased with a wetting rate +3.8%/day for 5 days. 

From the middle of December until the end of the monitoring 
period, no snow gliding was observed, although a significant in¬ 
crease of snow/soil interface temperature and moisture was re¬ 
corded at the end of March (wetting rate of +1.5% day for 12 
days). 
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Fig. 3 Snow depth and air temperature measured by the automatic snow and weather station (a); 77 (b) Mi (c) and snow gliding (d) in each plot 

recorded in both winters (2009-2010 and 2010-2011) when the study site was snow-covered. 


Snow forces 

In the Larch plot , the snow forces on stems with different 
diameters, calculated by Equation 3, as in Margreth et al. (2007), 
ranged from 3 to 290 N (Fig. 4a). Two maximum occurred at the 
middle of January ( H s = 60 cm and p snow =184 kg-m" 3 ) and at 
the beginning of March ( H s = 47 cm and p snow = 362 kg-m" ). 

We also calculated the snow forces as in Me Clung and Larsen 
(1989), corresponding to the single snow gliding episodes 
(beginning of December), when a movement of 3.9 mm/day was 
registered by the snow shoes. The necessary input values for the 
snow depth and density were taken from the snow pit dug on 9 th 
of December: H s = 25 cm and p snow= 250 kg-m' 3 . The resulting 
forces ranged from 11 to 284 N on the different stem diameters 
(Table 2). 

In the Spruce plot the snow forces on stems with different 
diameters, calculated by Equation 3, as in Margreth et al. (2007), 
ranged from 2 to 361 N (Fig. 4b). A maximum was observed at 
the middle of February, when H s was equal to 66 cm and p snow 
231 kg-m" 3 . 

During winter 2009-2010, it was not possible to calculate the 
snow forces as in Me Clung and Larsen (1989), as no snow 
gliding event occurred under Spruce. 
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Table 2. Snow forces values calculated according to Me Clung and 
Larsen (1989) for the different stem diameter classes under Larch (L) 
and Spruce (S) in winters 2009-2010 and 2010-2011. 


Date and snow 


Forces at different diameter classes (N) 


gliding rates 

0.02 m 

0.05 m 

0.1 m 

0.2 m 0.3 m 

0.4 m 

0.5 m 

2009-2010 a 




Larch plot 



December 4 th , 

3.9 mm/day 

11.37 

28.43 

56.86 

113.72 170.58 

227.47 

284.30 

2010-2011 




Larch plot 



November 17 th , 
1.4 mm/day 

7.38 

18.46 

36.92 

73.83 110.75 

147.66 

184.54 

January 20 th , 

0.4 mm/day 

2.55 

6.39 

12.77 

25.55 38.32 

51.10 

63.87 

March 03 rd , 

0.1 mm/day 

1.82 

4.54 

9.08 

18.17 27.25 

36.34 

45.42 





Spruce plot 



November 19 th , 
0.4 mm/day 

2.62 

6.56 

13.11 

26.22 39.33 

52.45 

65.56 


a In 2009-2010 no snow gliding occurred in the Spruce Plot. 
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Fig. 4 Snow forces calculated according to Margreth et al. (2007) at 
different stem diameters during winter 2009-2010 in the Larch (a) 
and Spruce plots (b). 

Winter 2010-2011 

Meteorological pattern 

In winter 2010-2011 the mean air temperature was -3.0°C with a 
minimum value of -14.3°C recorded on December 18th (Fig. 3). 
Winter 2010-2011 was characterized by frequent and longer 
mild periods than winter 2009-2010. The mean daily air 
temperature became constantly positive from March 22nd. 
Concerning precipitations, the automatic weather station 
recorded a cumulative precipitation equal to 298 mm of snow 
water equivalent (SWE). The maximum snow depth (156 cm) 
was reached on March 17th. 

The cumulative snowfalls were lower than the 75 years 
historical data (488 mm of SWE) reported by Mercalli et al. 
(2003): in particular, the difference was remarkable especially 
during January, February and April. 

Snow physical characteristics 

In the Larch plot , a maximum snow depth equal to 28 cm was 
observed on December 2nd. The snow density ranged between 
146 kgm" (on December 2nd) and 340 kg-nf (January 18th). 
The first snow profile showed a kinetic metamorphism with the 
presence of faceted crystals (FCso), which we found also on 
January 18 th . Moreover the dominant grain type in the basal layer 
was rounded polycrystals (MFpc) suggesting a previous presence 
of liquid water. On March 23rd a prevalence of cluster rounded 
crystals (MFcl) was observed. 

In the Spruce plot , the maximum snow depth was observed in 
the last snow profile (28 cm) on March 23rd. The snow density 
ranged between 230 kg-m" (on December 2nd) and 320 kg-nf 
(January 18th). The first snow profile revealed the beginning of 
kinetic metamorphism, with faceted crystals (FCso) as prevalent 
crystals; moreover a thin ice layer (IFil) was observed at the bot¬ 


tom of the snowpack. On March 23rd a prevalence of rounded 
crystals (RGlr) was observed. 

Temperature and moisture at the snow/soil interface 
The temperature sensors in the Larch plot were repeatedly dam¬ 
aged by rodents, therefore the measurements were interrupted. 
The recorded 7) values ranged only from November to December 
and from the end of January until the middle of February (Fig. 
3b). 

The moisture at the snow/soil interface M z (Fig. 3c) was char¬ 
acterized by 4 sharp increases: the first (34%) occurred from No¬ 
vember 14th to 16th (wetting rate: +12%/day); the second (34%), 
occurred from December 5th to 9th (wetting rate: +9%/day); the 
third (37%) occurred from January 13th to 16th (wetting rate: 
+5%/day); the fourth (32%), occurred from March 13th to 20th 
(wetting rate: +4%/day). On the contrary, M z strongly decreased 
(-25%) from November 21th until December 4th, with a water 
depletion rate equal to -1.8%/day. M z was positively correlated 
with the air temperature (r = +0.256; p < 0.01). 

In the Spruce plot the temperature at the snow/soil interface f 
(Fig. 3b) was constantly negative from November 24th until 
March 3rd, except three periods when it reached 0 °C (December 
8th, January 14th to 18th and February 6th to 11th). Two mini¬ 
mum values were reached, on December 18th (-4.6°C) and Janu¬ 
ary 21th (-4°C). When the snow/soil interface temperature in¬ 
creased, also the moisture M z did (r = +0.650 p <0.01), and the 
maximum values, 42%, 49% and 33%, were reached on Novem¬ 
ber 16th, January 15th and February 9th, respectively (Fig. 3c). 
Instead, the minimum water content measured in the Spruce plot 
(0%) was reached at the beginning of December when the 
snow/soil interface was frozen. 

Snow gliding 

In the Larch plot , several snow gliding events were measured 
with a cumulative movement of 77 mm (Fig. 3d). The first snow 
gliding event was recorded after the first snowfall (25 cm) at the 
end of November. From November 13th to 23rd a cumulative 
movement of 43 mm was recorded. In this period, Ti was posi¬ 
tive (daily average 7)= 1.4 °C) and a remarkable snow/soil mois¬ 
ture increase (+ 7.4%/day for 3 days) occurred. The highest snow 
gliding rate (1.4 mm/day) was measured on November 17th. In 
the snow profile dug closer in time to this snow gliding event 
(December 2nd) the dominant grain type in the basal layer was 
rounded polycrystals (MFpc), indicating that melt-freeze meta¬ 
morphism occurred, with the presence of some liquid water. Un¬ 
til the middle of January, the snow gliding rates decreased, with 
a cumulative value equal to 5.8 mm. From January 19th till April 
4th other snow gliding events were observed (with a cumulative 
movement of 29 mm): the highest daily movement was recorded 
on January 20th (0.4 mm/day) and on March 3rd (0.1 mm/day), 
after a significant increase of M z (respectively +2.4/day for 7 
days and +8.7 in one day). Especially on January 18th an in¬ 
crease in snow moisture across the snow pack was observed: the 
basal layer contained a significant amount of liquid water. No 
snow gliding was observed from February 6th until 16th because 
the study plot became snow free. 
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In the Spruce plot , only a single snow gliding event was ob¬ 
served, after the first snowfall, with a cumulative displacement 
of 16 mm, from the 13th to the 25th November. In the snow pro¬ 
file dug closer in time to this snow gliding event (December 2nd) 
a thin ice layer (IFil) was observed at the bottom of the snow- 
pack indicating a previous presence of some liquid water at the 
snow/soil interface. The highest snow gliding rate (0.4 mm/day), 
was measured on November 19th, after a significant increase of 
the volumetric moisture (+4.75/day for 3 days). From the end of 
November until the end of the monitoring period, no snow glid¬ 
ing was observed, despite several increments at the snow/soil 
interface moisture were recorded. 

Snow forces 

In the Larch plot , the snow forces on stems with different 
diameters calculated by Equation 3 (as in Margreth et al., 2007) 
ranged from 0.5 to 29 N (Fig. 5a). The snow forces decreased 
from the beginning of December ( H s = 28 cm; p snow= 146 
kg-m" 3 ) to the beginning of January ( H s =10 cm; p snow = 340 
kg-m" ) and then remained almost constant until the end of the 
season. Also according to Me Clung and Larsen (1989), the 
snow forces decreased across the season since the maximum 
snow gliding rate, 1.4 mm/day, was observed on November 17th 
and the minimum one, 0.1 mm/day, was measured on March 3rd. 
The resulting forces ranged from 2 to 184 N (Table 2). 

In the Spruce plot , the snow forces on stems with different 
diameters calculated by Equation 3 (as in Margreth et al., 2007) 
ranged from 0.6 to 21 N (Fig. 5b). The snow forces were almost 
constant until the beginning of January ( H s =11 cm; p snow = 
320 kg-m" ) and then started to increase until the end of March 
(H s = 28 cm; p snow 260 kg-m' 3 ). We also calculated the snow 
forces as in Me Clung and Larsen (1989), in correspondence of 
the single snow gliding episode (second half of November), 
when a movement of 0.4 mm was registered by the snow shoes. 
The necessary input values for the snow depth and density were 
taken from the snow pit dug on December 2nd: H s = 20 cm and 
p snow= 232 kg-m" . The resulting forces ranged from 3 to 66 N 
on the different stem diameters (Table 2). 

Comparison between Larch and Spruce plots 

Snow physical characteristics 

Tree species affect only some physical characteristics of the 
snowpack. In fact, the crystal type and dimension were similar 
both under Larch and Spruce stands. Instead, the snow depth 
evolution had a different pattern for the two sites: in mid winter, 
snowfalls generated a deeper snowpack in the Larch plot, 
because of the reduction of snow interception due to the absence 
of leaves in the winter period; on the contrary, in spring, the 
shadow effect due to the evergreen cover maintained a deeper 
snowpack and increased the snow duration under the Spruce plot. 
The same effect of crown interception was found by Gubler and 
Rychetnik (1991) and Motta (1995). 

A higher average snow density was observed under Spruce, 
especially after snowfalls with mild air temperature, maybe due 
to the liquid water input from the melted snow intercepted by the 
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tree crowns. The snow density under Larch gradually increased 
and finally exceeded the values recorded under Spruce in spring, 
because of a more pronounced effect of solar radiation. 
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Fig. 5 Snow forces calculated according to Margreth et al. (2007) at 
different stem diameters during winter 2010-2011 Larch (a) and 
Spruce (b) plots. 


Temperature and moisture at the snow/soil interface 
The mean snow/soil interface temperature was higher in the 
Larch than in the Spruce plots, especially during fall and spring. 
An explanation could be that the greater snow accumulation un¬ 
der Larch during early winter, when the soil was unfrozen, could 
have reduced the soil cooling rate, while the higher snow depth 
in the Spruce plot in late winter and spring could have deter¬ 
mined a delay in the soil warming. Several authors (e.g. Brooks 
and Williams 1999; Freppaz et al. 2008) found that a snow depth 
approximately equal to 30-40 cm accumulating early in the win¬ 
ter season prevents soil from freezing. 

Concerning the moisture at the snow/soil interface, we regis¬ 
tered larger fluctuations in the Larch than in the Spruce plot. The 
reason might be that a lower crown cover allowed the radiation 
to reach the snowpack surface. This energy input increased the 
snow water content at the surface, which, due to a thin snowpack, 
was able to percolate downwards to the ground. On the contrary, 
in the Spruce plot, the larger crown cover lead to a lower influ¬ 
ence of solar radiation and probably prevented the wetting front 
from reaching the soil. 


Snow gliding 

In both years the total snow gliding was higher under Larch than 
under Spruce. This fact was probably due to a lower snowpack 
basal friction under Larch than under Spruce, caused by great 
fluctuations of moisture at the snow/soil interface. The cumula- 
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tive snow gliding did not exceed 100 mm, as found also by 
Holler (2001) for a dense forest. Our results support the hypothe¬ 
sis that snow gliding is generally low when the canopy density is 
relatively high. 

Snow forces 

Generally, lower forces were calculated below Spruce than be¬ 
low Larch, due to the larger snow accumulation recorded in the 
latter site. We did not find any damages to young trees in both 
plots. 

General discussion 

From our results, the snow/soil interface moisture significantly 
influenced the snow gliding processes. In particular, soil mois¬ 
ture at the snow/soil interface was related to the snow/soil inter¬ 
face temperature and both parameters (Mi and Tj) were strongly 
regulated by air temperature. It seemed that the wetting rate at 
snow/soil interface was strongly influenced by air temperature 
and snow depth with a contrasting effect. In particular, we found 
that during early winter, when the snow cover was thin and ho¬ 
mogeneous (mainly characterized by fresh snow) the volumetric 
moisture at the bottom increased in less than 24 hours under both 
plots after an increase of air temperature, as found also by Jones 
et al. (2001). Instead, in mid winter, when the snowpack was 
thicker and less homogeneous, after an increase of air tempera¬ 
ture, the snow/soil interface moisture increased because of either 
a wetting front from the surface (which reached the bottom of the 
snowpack within a maximum of 4 days) and an increase of the 
snow/soil interface temperature (which reached 0 °C). 

The volumetric moisture increase of a relatively dry snow/soil 
interface (10%-20%), seemed to be the most important factor 
influencing the snow gliding rates: in particular snow gliding 
was observed with wetting rates higher than 1.2%/day. The in¬ 
crease of volumetric moisture at the snow/soil interface might 
reduce the snowpack basal friction as suggested by several au¬ 
thors (In der Gand 1966; Clarke and Me Clung 1999; Holler 
2001; Margreth et al. 2007). The snow gliding events were 
mostly recorded, under both plots, during early winter and after 
the first snowfalls, when the soil was unfrozen, as observed also 
in long term studies, though conducted in open field (In der Gand 
and Zupancic 1966; Hoeller 2001). The first snow reaching the 
ground rapidly melted on a relatively warm soil (7) >0.5 °C), 
thus determining an increase of moisture at the snow/soil inter¬ 
face, where the dominant crystal type was slush. During winter 
2010-2011, snow gliding events were observed also during 
spring time. These phenomena were related to a sudden sharp air 
temperature increase, which consequently increased the snow¬ 
pack moisture, followed by a strong air temperature decrease, as 
observed by McClung (1999). We could hypothesize that the 
snowpack became stiff and glided, as a block, on the liquid water 
film previously formed at the soil surface. 

The snow forces calculated in this study were lower than the 
values reported in other studies conducted in the Alps (Holler et 
al. 2009). The differences could be attributed to the lower coeffi¬ 
cients we used, because our study site was located under a ma¬ 


ture forest stand, while the study site considered for example by 
Holler et al. (2009) was characterized by young trees. Moreover, 
our study site was characterized by lower snow density and 
thickness than the site studied by Holler et al. (2009). 

With our data, the application of the method by Margreth et al. 
(2007), originally elaborated for the dimensioning of the cable- 
way masts, gave lower values than the method by Me Clung and 
Larsen (1989), and the differences between the two methods be¬ 
came more significant with larger stem diameters. This differ¬ 
ence highlights the need to directly measure the snow forces, by 
e.g. using load cells on the tree trunks, in order to validate the 
results obtainable with the two different methods. 

The fact that we did not detect any damages to young trees, 
such as seasonal uproots and branch breakages, is in agreements 
with the low snow forces obtained with the two different meth¬ 
ods. Our calculated snow forces were lower than the value of 950 
N, found by Holler (2009), through pulling test, as the minimum 
force necessary to cause uprooting of young coniferous trees 
(Larix decidua, Picea abies e Pinus cembra ), with stem diame¬ 
ters ranging from 2.5 to 5 cm. 

However, the snow characteristics (e.g. density) and meteoro¬ 
logical conditions (e.g. sharp increase of air temperature) re¬ 
corded during the two winters considered in this study, caused 
significant snow gliding events, mainly during early winter, 
which contributed to determine non negligible snow forces on 
trees, though not producing damages. 

Since the Climate change is expected to determine milder win¬ 
ter and more unstable snow cover in temperate mountain regions 
(Cooley 1990; IPCC 2007), it is possible to speculate a change 
also in the snow gliding processes. In particular, these areas 
could face a decline in snow cover (late snowpack accumulation 
and advance of snowmelt timing), an increase of snow density 
and more frequent episodes of rain on snow (Leung et al. 2004; 
Ye et al. 2008; Casson et al. 2010). Hence, the basal friction of 
the snowpack could be reduced by the liquid water on the basal 
portion of the snowpack and consequently amplify the rate of 
snow gliding. This fact might generate in future higher pressure 
on stems diminishing trees stability and therefore the forest pro¬ 
tective function. 

Conclusion 

In this work we evaluated the relationship between snow physi¬ 
cal properties, microclimatic conditions at the snow/soil interface 
and the process of snow gliding during two winter seasons under 
two distinct forest covers, by means of continuous measurements 
and detailed field observations. 

In particular we observed that: 

(a) The snowpack physical evolution (e.g. snow density, snow 
depth) during the whole winter season was influenced by the tree 
species, due to the different crown interception of snowfall, and 
the different degree of exposure of the snowpack to the incoming 
solar radiation; 

(b) The wetting rate at the snow/soil interface, associated to air 
temperature fluctuations, seemed to be the most important factor 

^ Springer 




642 


Journal of Forestry Research (2013) 24(4): 633-642 


influencing the snow gliding rates. The snow gliding rates were 
higher in early winter, when the first snowfall covers an unfrozen, 
warm soil. In this period, snow gliding was observed after the 
wetting of a relatively dry snow/soil interface; 

(c) The highest snow gliding was measured under the Larch 
stand probably due to a lower basal friction at the snow/soil in¬ 
terface; 

(d) The calculated snow forces were very low under both for¬ 
est stands, if compared to other studies, mainly because the snow 
depth and the snow density were scanty during the study period. 
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